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Jet Function with a Jet Algorithm in SCET 
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The jet function for the factorized cross section e + e~ into dijets is given as a function of the jet 
invariant mass s and with a generic jet algorithm at 0(a 3 ). We demonstrate the results using the 
Sterman- Weinberg algorithm and show that the jet function is independent of the energy fraction 
/3 of the soft radiation. The anomalous dimension has the same form with and without the cone 
half-angle 8. The dependence of the finite part of the jet function on the cone angle is given. 



I. INTRODUCTION 

Hadronic jets feature in many final states of interest in modern collider experiments. They form a significant 
Standard Model background for many new physics processes and also provide probes for QCD interactions at several 
different scales. In order to disentangle the effects of these different momentum scales, factorization is needed. Factor- 
ization theorems make it possible to separate process-dependent perturbative physics from universal nonperturbative 
effects. Soft-collinear effective theory (SCET) [3-lE| provides a framework for deriving factorization theorems while 
systematically resumming large logarithms to all orders in perturbation theory and including power corrections to 
any desired accuracy. The power expansion is performed in terms of a parameter A, which characterizes the ratio of 
the transverse and the collinear momentum in a jet. 

The cross section for e + e~ into dijets can be factorized schematically as follows 

— 4^ =H{J n ® Jn®S), (1) 

do ds as 

where s and s are the invariant masses of the jets, do is the Born cross section, H is the hard coefficient obtained by 
matching SCET to QCD, J„/ fi are the jet functions in n and n directions, S is the soft function, and the convolution 
is in the small light-cone momentum component for each jet. We study a cross section differential in s and s because 
after specifying the total energy and the direction of a jet, the invariant mass is the next natural observable to consider 
in order to find out more about the structure of the jet. Factorization theorems for e + e~ colliders have been derived 
in SCET for event shape variables using a hemisphere jet algorithm For hadron colliders, factorization theorems 

have been studied at a general level [8| as well as for threshold resummation and for isolated Drell-Yan [l8[ 

defined as having no central jets. 

To match experimental results more precisely, a theory calculation should define jets using the same jet algorithm 
that is used in a given experimental analysis. The present work is a first step towards the final goal of proving a 
factorization theorem for the dijet cross section given above with a realistic jet algorithm. Instead of attacking the 
problem all at once, we derive a formalism to calculate the complete jet function to 0(a s ) and to leading order in power 
counting. Our results apply for any jet algorithm that can be formulated in terms of theta functions depending on 
the momenta of the final state particles and that can be regularized using dimensional regularization. As an example, 
we look at the Sterman- Weinberg algorithm (SW) which has been considered using SCET in Ref. [20-23]. We 
calculate the full jet function J„, including the finite part, as a function of the invariant mass of the jet. After an 
expression for the soft function has been calculated, our result can be used to derive a factorized cross section e + e~ 
into dijets with the jets defined using the SW algorithm. 

The SW algorithm defines a dijet event as one where all but a fraction f3 of the total energy is contained within a 
pair of oppositely directed cones of half-angle S The definition can also be extended for more jets if one specifies 
a way to determine the direction of the different jet axes. Very recent work by Ellis et al. [24| shows that in the case 
of more than two jets, consistent factorization requires the jets to be well separated and the radiation outside the jets 
to be soft. Hence, we expect a factorization proof for N jets defined with SW algorithm to take the schematic form 



1 dg bW (<5,/3) 
ctq dsi . . . dsN 



N 



= h n j » sw ^ sw , (2) 



with 0(6,(3) power corrections. For the expected power corrections to be subleading in power counting, it is natural 
to take S,/3 ~ A. As will be seen later, to leading order in (3 the final result for the jet function is independent of (3, 
even if we pick S ~ A, (3 ~ A 2 . It would seem logical to also consider 5 ~ A 2 , (3 ~ A but as has been pointed out in 
Ref. [IH, one should require sm(S) > f3/(l — (3) in order to preserve the back-to-back orientation of the two jets. To 
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lowest order in the parameters, this corresponds to S > (3. Finally, one could consider 5, (3 ~ A 2 but then it would 
be more natural to define a new expansion parameter A' = A 2 and set S,/3 ~ A'. Thus, the natural possibilities to 
consider are <5^/3^Aor<5^A, (3 ^ X 2 and our results for the jet function are valid for both cases. For the sake of 
definiteness, we will take S ~ (3 ~ A with the understanding that sin(<$) > (3/(1 — (3). 

Ref. [24[ also provides a useful cross check for our results. They have demonstrated the consistency of factorization 
of jet observables in exclusive multijet cross sections for both cone and cluster type algorithms and given expressions 
for the anomalous dimensions of the hard, jet, and soft functions. The present work agrees with their result for the 
anomalous dimension of the jet function. We also derive the finite contribution, which has not been given in the 
literature before. 1 

In Sec. |n] we review the necessary elements of SCET and present our results for a generic jet algorithm. In Sec. IIHI 
we express the SW jet definition in terms of the components of the momenta of the final state particles. In Sec. IIVI 
we present the main results of our calculation. We plot the renormalized jet function both as a function of s and 
integrated over s. For comparison, the jet function is also plotted without a jet algorithm. We conclude in Sec. fVl 



II. JET FUNCTION WITH GENERIC ALGORITHM 



Our goal is to derive a formalism to accommodate any jet algorithm that can be expressed in terms of phase space 
cuts. In this section we outline the derivation for an expression for the jet function J rl (s). We start by reviewing 
some conventions for notation. The direction of the jet axis is denoted by a unit vector ft. We choose a coordinate 
system such that n — (0,0,1). It is convenient to work in the light cone coordinates with basis vectors = (l,n) 
and — (1, — n) satisfying n 2 = n 2 = 0, n-n — 2. Then any momentum can be decomposed as 

P 1 * = n-p — + n-p — +p' J L = (n-p, n-p, p ± ) = (p + , p~, p ± ) . (3) 

The momentum of a collinear parton i in n-direction scales as pi ~ Q(A 2 , 1, A) so we can write p®, p^ = + 0(A). 
The center-of-mass energy is denoted by Q and the total jet momentum is p n + r n = (Q,r + ,0), where p n is the 
large label momentum and r n is the smaller residual momentum. Thus, the invariant mass of the jet is given by 
s = (p n + r„) 2 — Qr+ ~ Q 2 A 2 . Because of the specific observable we consider, namely the jet invariant mass, we are 
free to use the coordinate system specified above, which removes all dependence on the total transverse momentum 
Pn ■ For a general observable, we would also need to include p^. To define a jet algorithm, we also need to consider 
the momenta of the final state partons. At 0(a s ), momentum conservation and the on-shell condition for the final 
state partons make it possible to write all the momenta in the problem in terms of the jet momentum components 
Pn = Q and together with the gluon four- momentum p g . 

The jet function can be written in terms of the gauge invariant quark jet field x„ = £n, and it is also convenient 
to use Xn,uj = S(lo — n-V)(W^ n ), where n-V gives the large label momentum of the combination M^J^. The collinear 
Wilson line is defined as 



W n (x) = ^2 ex P 



perms 



n-V 



n-A n (x) 



(4) 



where A n is the collinear gluon field. 

We will use the term inclusive jet function to refer to the case where no jet algorithm is applied and denote it by 
J n 1Tlc ^ . It is discussed in Ref. 0, [53] and we define it as 

J " inC)(SlAi) = 8^Q DiSC J ddxelrn ' Xil (°\ T ^, Q mXn(x)}\0) 

1 

= 8^ 



^£ Jd d xe^tT(Q\ii Xn {x)\X n )(X n \x^ Q (m), (5) 



where N c is the number of colors, d = 4 — 2e, the trace is over color and spin, and T stands for time ordering. Because 
of charge conjugation symmetry, the antiquark jet function J n does not have to be considered separately. We write J n 



Right after this paper appeared, a second paper by Ellis et al. [25j came out. It is discussed in a note added and an Appendix at the 
end. 
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FIG. 1: The Wilson line diagrams (a) and (b) give identical contributions. Diagram (c) vanishes in Feynman gauge. Diagram 
(d) gives the wavefunction renormalization contribution. Diagrams (e) and (f) show the real cuts and diagrams (g) and (h) 
show the virtual cuts. The mirror image of (h) gives the other virtual cut of (d). 



as a sum over final states as in Ref. @ in order to implement a jet algorithm. The final states are restricted according 
to a constraint function F(ai), which defines an algorithm in terms of parameters <Zj and depends on the momenta 
of the particles in \X n ). Thus, F(ai) is also a function of the operators pj which have the final state momenta pj 
as eigenvalues. In this section, we will work with a generic F{ai 1 f>j) and in Sec. IIIII we will specialize to the SW 
algorithm. Inserting the constraint function gives the algorithm-dependent jet function 

J " (s,Ai) = 8^Q^ /d^e^-tr^l^xn^)^^,^)!^)^!^^^)^). (6) 

Eq. ([6]) is valid at any order in a s but the functional form of F(cii,pj) changes from order to order. Following the 
discussion of Hornig et al. on angularity jet functions in Ref. Q , we implement the needed phase space restrictions 
on the final states by introducing an ^-discontinuity, where the standard Cutkosky cutting rules are modified by 
inserting a factor of F(a,i,pj) into the cut propagators. The diagrams that contribute to J„ at one loop are shown in 
Fig. [T] (a) - (d). We note that diagram (c) vanishes in Feynman gauge. At this order, we can cut through the loop or 
a single quark propagator, which we call "real" and "virtual" cuts, respectively. Diagrams (e) and (f) give the real 
cuts and (g) and (h) show the virtual cuts. The virtual cuts are independent of the algorithm and are contained in 
a proper interpretation of the inclusive jet function. Hence, the algorithm-dependent contribution is given by taking 
the real cuts and inserting F(cii,pj) to the cut propagators. In order to use the known results from the literature 
more conveniently, we add and subtract the inclusive jet function on the second line. 

J " (S) = 8^Q D ' 1SCF /d^^»-tr<0|T{x„ )Q (0)^Xn(^)}|0> 

= 4 inc) (*) + ^g Disc /d'W— tr(0|T{x„,g(0) {F{a u p 3 ) - 1) fe n (x)}|0) = J^OO + Ajf, (7) 

where we have introduced notation A j£ for the algorithm dependent contribution. Note that setting F = 1 in Eq. ([7]) 
leads to AJ,f = leaving only the inclusive contribution, as expected. To find an expression for AJ„ , we use the 
modified cutting rules to add a factor of F(ai,pj) and to replace the gluon and quark propagators in the loop by delta 

functions S (p-pj + (pj) 2 ) and 8 (p-pj + (pj) 2 ) . Using the momentum conservation relation p q + p g = (Q, r+ 0) 

and the above delta functions, the phase space integrals over the gluon momentum components pj and pj can be 
performed to give the relations 

Pg p - > Wg) Q 2 ■ W 

These relations must be used when writing the constraint function. Finally, the calculations can be simplified by a 
change of variables y =p~/Q, 

In a factorization theorem the contribution of soft quarks and gluons is encoded in the soft function. In order 
to avoid double-counting when the loop momentum of a collinear field in the jet function becomes soft, a zero-bin 
subtraction must be performed [28(. The naive collinear result without the subtraction is obtained by summing over 
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the contributions from the real cut Feynman diagrams in Fig. Q] diagrams (e) and (f ) , (and counting diagram (e) 
twice to account for its mirror image). This gives 

A J,f (.) = (£) ^f^y 1(1 - „)- (H-^l + y(d - 2)) (F(a,i,y) - 1), (9) 

where Eq. (jHJ has been used to write f 1 as a function of y and the algorithm parameters a«, the tilde denotes that 
the zero-bin subtraction has not been performed, and A(e) = 1 — ^rj — h C(e 3 )- For the SW algorithm, we will see 
that after the zero-bin has been subtracted, AJ„ is finite as e — > so we will eventually take A(e) = 1. There is a 
zero-bin contribution both for the gluon and the quark becoming soft but only the former contributes at leading order 
in power counting. Furthermore, as the soft gluon wavefunction renormalization vanishes in Feynman gauge, we only 
need to consider a zero-bin for the gluon in Fig. [T] (e). To obtain the zero-bin result, we assign scaling p g ~ QX 2 [28| 
to all components of the gluon momentum. The effect of the zero-bin scaling to algorithm constraints is discussed in 
detail in Appendix [A] The end result is that instead of Eq. (j8]), the replacement rules for the zero-bin piece are 

pt = %> (pi) 2 = -PgPt = -^-> p*=Q> pJ = °. pt = °- ( 10 ) 

Using these relations, the zero-bin contribution to A Jjf is 

A^foOO = ^T A (^ fay ^(Fo(a u y) - 1), (11) 

where Fo(a,i,y) denotes that the constraints are written using Eq. (|10[) . For a very inclusive jet algorithm such 
as the hemisphere algorithm, the zero-bin contribution remains a scaleless integral that vanishes in dimensional 
regularization. However, a more restrictive algorithm can introduce a scale that gives a nontrivial zero-bin subtraction 
and we will see that this is what happens for the SW case. 



III. STERMAN- WEINBERG ALGORITHM 



As explained in the introduction, the SW algorithm defines a dijet event in terms of the cone half-angle 6 and 
energy fraction of the soft radiation (3. We will keep only the leading order in 5 ~ j3 ~ A, including both power and 
logarithmic dependence. When we apply the SW algorithm to the jet function, the question at one loop is whether 
the n-collinear quark and gluon create one or two jets. If both partons lie within 5 from the jet axis or if one of 
the partons has energy Et < (3Q, only a single jet is produced. As stated before, a soft quark would contribute 
only at higher order in the power counting so we have two regions of phase space to consider: a "cone" region where 
both particles lie inside the cone and an "outside-cone" region where the gluon is not inside the cone and has energy 
Eg < f3 Q. We take the jet axis to lie along the n-direction and then the cone region corresponds to the momenta of 
both partons lying within 5 of n. 

The invariant mass s of the jet is related to how far apart the final state particles are spread. At one loop and 
using the on-shell condition for massless partons, we find that s = PgPq sin 2 (</>/2) ss p~p~(cf>/2) 2 where <fi is the angle 
between the quark and the gluon. If no jet algorithm is applied, in the center-of-mass frame momentum conservation 
implies s < (Q/2) 2 . In the SW algorithm with the gluon inside the cone, the maximum value is reached when 
Pg = Pq = an d 4> — 2<5, which gives s < (QS/2) 2 . This constraint is satisfied by the naive contribution to the jet 
function which only has support for s < (QS/2) 2 . However, the standard approach to the gluon zero-bin subtraction is 
to utilize the result obtained in the zero-bin limit everywhere, and hence there is no upper limit on p~ or on s. It may 
be possible to use the freedom in defining a zero-bin subtraction scheme to define a non-minimal subtraction where 
one would obtain a jet function which is nonzero only for s < (Q5/2) 2 . We leave this question to future work. When 
the gluon is outside the cone, the maximum is reached when y = 2/3 and <f> = ir, which gives s < 2/3 Q 2 . However, 
because of the power counting s ~ QA 2 , this condition is always satisfied to leading order. 

The constraints for the SW algorithm are illustrated in Fig. [5] in terms of the variables p~ and . Fig. [5] (a) 
describes the phase space for the naive collinear contribution. If there were no jet algorithm, pj^ would be integrated 
all the way to infinity in the full range < p^ < Q. For the SW algorithm, the region defined by p~ < 2/3 Q and 
2\pg L \/p~ > S corresponds to a gluon being emitted outside the cone, and the triangular region at the bottom of the 
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figure gives the cone contribution. The corresponding constraint function -Fsw is given by 

V 




4s + Q 2 S' 2 



^-t&n8je(pQ-p° g ) 
y ) + 9(2(3-y) + 0(X,S,(3), (12) 



where Eq. ® was used in getting to the second line, and the first thcta function in the second term is always satisfied 
at leading order in the power counting. This is because As/ {As + Q 2 S 2 ) ~ A but y ~ A in the second term due to 
/3 ~ A. We note that the two theta functions in the first term imply 9(Q 2 S 2 /A — s) and thus limit the maximum 
allowed jet mass. 

The zero-bin phase space is shown in Fig. [2] (b). Without the algorithm, the integration region would extend to 
infinity for both p~ and p^ . As explained in Appendix [XJ the zero-bin scaling only affects the jet algorithm through 
the conditions in Eq. (|10p. This results in a zero-bin constraint function 

F , sw = 6 (tan* - +8^ - UnSj 9 (fiQ-fl 

= e {y--^l) + W-y) + 0(A,<5,/3), (13) 

where we again used the fact that As/(Q 2 8 2 ) ~ A to eliminate the first theta function in the second term. (This result 
for the constraint function in the zero-bin region agrees with Ref. [25s].) We note that the conditions from Eq. (fT0|) 
eliminate the second theta function in the first term. This theta function constrained the quark to be inside the cone. 
Its absence can be understood physically by remembering that in the gluon zero-bin, the quark carries all of the label 
momentum and is automatically inside the cone. In this case, there is no upper limit on the jet mass s. 



IV. RESULTS 



A. Algorithm-dependent Contribution 



We are now ready to apply our general result for the case of the SW algorithm. For both y and s/ pi 2 , we have used 
the following distribution identity 



o(y) <%) 



y 



l+e 



%) 



V 



9(y)Hy) 



0(e 2 



(14) 



where [9(y)\n n (y)/y] + denote plus functions, which we define such that they give zero when integrated from to 1. 
The different possible definitions and their relationship with one another are discussed for example in Appendix B of 
Ref. [13. 

Combining Eqs. ([9]) and (|12[) we find for the naive collinear contribution 



Aj sw (s) =^L A (e)9 (>-i±)^ ( s, 



4--(hi2/3 + iln. , 
e 2 e| 2 V Ajj? 



2l„-U-hr (9£ 



^ 2 9(s) 



In 2/3 + In 



12x2 



Afi 2 



6 s 



4s + Q 2 S 2 



p 2 9{s)\n(s/ix 2 



An 



4s ^ ± 
Q 2 ' 1 



pi 2 9{s) 



- +ln2/3 
4 



(15) 
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FIG. 2: Phase space regions for the SW algorithm. Naive collinear (a) and zero-bin (b). The region at the left-hand side of the 
figures corresponds to a gluon outside the cone with energy E g < /3Q. The region at the bottom of the figures describes the 
region where both partons are inside the cone. For the zero-bin contribution, pj is integrated all the way to infinity. 



Eqs. and (|13|) give the corresponding zero-bin contribution 



AJ,f„ w (») 




In 2/3 



+ 



a s C F ( 

47T 



4s 

Q 2 









s 



1, (Q 2 5 2 



/i 2 6>(s) In (s/m 2 ) 
s 

2x2 



In 2/3 + In 



Q S 
<±H 2 



/z 2 0(s)ln (s/^ 2 ) 



The total algorithm-dependent part A J, 

a s CF 



AJ r f w ( s) 



4tt 



- ^7 



is given by the difference of Eqs. (| 15[) and (TTB1) and reads 

CV S C_F 



4s 



24 



4s + Q 2 <5 2 4tt 



ii-<5 2 

Q 2 



3 4, 
- + - In 

s s 



4s 
^J 2 



(16) 



(17) 



Using the definition of the plus function [30], it can be seen that when [Q(y)\i\ n {y)/y\ + is integrated against a 
function /(y) such that /(0) = 0, we can ignore the plus function prescription because it only makes a difference at 
y = 0. Hence, we do not need to use plus functions in Eq. (fT?) . We note that the outside-cone region for the gluon 
cancels between the collinear and the zero-bin contributions and hence there is no /3-dependence in the final result. 
This cancellation has been discussed in the literature [23| and has to take place since the purpose of the zero-bin 
subtraction is to remove any soft contribution from the collinear diagrams and the gluon must be soft to be outside 
the cone. We have also grouped the terms according to whether s is limited to be less or greater than (Q5/2) 2 . 



B. Full Sterman- Weinberg Jet Function And Anomalous Dimension 



We will now give an explicit expression for the full SW jet function. To do so, we need the inclusive jet function 
J^ 1110 - 1 , which can be found in the literature (2(| H?} • Using the notation of Ref. @ , it can be written as 



4 inc W) = s(s) + 



a s C F 

47T 



6(8) 



7-n z 



)-- 


u 2 e( s y 




-3 


4 

+ — 


> 2 6 


(s)ln 


) V 2 


s 






M 


s 



(18) 



The SW jet function is given by the sum of Eqs. (fT?)) and (|18p . Since the jet algorithm does not modify the 1/e-poles, 
the anomalous dimension is 



7i w (^) 



{- 


u 2 9( s y 


-6 5(s)| 




s 





(19) 
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0.00 0.01 0.02 0.03 0.04 0.00 0.01 0.02 0.03 0.04 

S [2 2 /4] S [2 2 /4] 

FIG. 3: The renormalized jet function in units of a s (fJ,)CF /kQ 2 , where we choose /i = 0.1 x Q/2. The SW (solid black line) 
and the inclusive (dashed purple line) jet functions are plotted as well as the difference of the two (dash-dotted blue line) for 
S = 0.2 (left) and 5 = 0.15 (right). 



The renormalized SW jet function in MS-scheme is given by 



J™ n ( s , M ) = S(s) + ^\S( S ) (7 - n 2 ) - A 



4tt 
4s 



24 



Q 2 J 4s + Q 2 S 2 



V0(s) 



4s 

Q 2 



H 2 0(s)ln(s/iJ 2 ) 



3 4 / 4s 

s + 7 ^Q^ 2 



(20) 



The most important result of this work is the renormalized SW jet function in Eq. (f2"0")) . It is plotted in Fig. [3] for 
S = 0.2 and 5 = 0.15 while keeping s ^ 0. We have chosen the jet scale to be /i = 0.1 x (Q/2). The SW and the 
inclusive jet functions are denoted by the solid black line and the dashed purple line, respectively. The jet function is 
continuous, but the derivative of the algorithm dependent contribution given in Eq. (|17p changes sign at s = (Q5/2) 2 , 
where one theta function turns off and the other turns on. We note that the algorithm contribution denoted by the 
dash-dotted blue line is always positive. 

Comparing the algorithm contributions in Fig. [3] for the two values of the cone angle, one can see a change in 
shape. For S — 0.15, small values of s — which correspond to narrow jets — contribute more than large values 
of s. For 5 = 0.2, the algorithm contribution is flatter, signifying a more equal contribution from jets of different 
size. This makes physical sense since a broader cone allows contributions from wider jets. Unlike the shape, the 
overall normalization of the algorithm contribution does not follow from physical intuition — we might have naively 
expected that increasing the cone angle would give a larger value for the magnitude of the jet function. However, 
care must be taken in assigning physical meaning to the normalization of a jet function computed in an unphysical 
subtraction scheme (MS), which is also illustrated by the fact that the inclusive jet-function has a zero at finite s. 
Different renormalization schemes include different constant pieces in the jet function and make its magnitude scheme 
dependent. Additionally, it should be remembered that a jet algorithm also affects the soft function, which must be 
combined with jet functions to see the full S dependence of the cross section. 



It is also interesting to integrate Eq. ([20)) over s up to s r 
and it is given by 



We call this the integrated jet function Jn 



SW(s 



V) 



°n,rcn Xt 1 ) 



= 1 + 9 S 2 - 



4s r 



Q 2 



ds J„ rcn (s,/z) 



as ° F h-n 2 -3ln /S ' 



4?r 



4s n 



<v 2 )^ 7-7r 2 -31n 



Q 2 S 2 
4/i 2 



2 ln A 



2 In 



2 5 2 



4/i 2 



6 In 



4s r 



Q 2 S 2 



2 <P 



61n2 + 41n| ^ ) In 1 4Smax 



Q 2 S 2 



(21) 



The renormalized NLO contribution to the integrated jet function is plotted in Fig. [4] as a function of s max and we 
have again chosen the jet scale to be fx = 0.1 x (Q/2). On the first row of Fig. 0] we plot the jet function for the 
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0.00 0.01 0.02 0.03 0.04 0.05 0.06 



Smax [0 2 /4] 

FIG. 4: The renormalized NLO contribution to the integrated jet function is plotted in units of a s (/i)CF/47r, where we choose 
/i = 0.1x Q/2. On the first row we plot the SW (solid black line) and the inclusive (dashed purple line) integrated jet functions 
as well as the difference of the two (dash-dotted blue line) for 8 — 0.2 (left) and S = 0.15 (right). On the second row we compare 
the integrated SW jet functions for the two values of S. 



SW case and the inclusive case, as well as the difference of the two, for S = 0.2 and 6 = 0.15. Since the difference 
between the SW and the inclusive jet functions was always positive in Fig. [3J we can see that the difference between 
the integrated SW and inclusive jet functions in Fig. [4] increases monotonically as a function of s max . On the second 
row of Fig. [5J the integrated SW jet functions for the two different values of 5 are compared. 



V. CONCLUSIONS 



Using a generic jet algorithm, we have derived an expression for the jet function j£ (s, /Lt) as a function of the jet 
invariant mass s, valid up to 0(a s ) and to leading order in the power counting. Expressions for the naive collinear 
and the zero-bin contributions have been given. We have demonstrated the general result by calculating the jet 
function together with its anomalous dimension for the Sterman- Weinberg algorithm, which is parameterized by the 
cone half- angle S and the energy fraction of the soft radiation ft. The anomalous dimension of the jet function is not 
modified by the jet algorithm. It was shown that all ft dependence is canceled because the gluon outside the cone must 
be soft and hence its contribution is removed by the nontrivial zero-bin subtraction. Our result for the anomalous 
dimension agrees with that reported in Ref. [24[. We have also calculated the finite part of the jet function. 

The renormalized jet function has been plotted as a function of the jet invariant mass to illustrate the difference 
between the SW and the inclusive case. We have also defined the integrated SW jet function Jjf w ( Smax ' by 
integrating J„ w (s,/z) up to s max . We have shown how J T f w (/z) changes as a function of S and how it differs from the 
inclusive case. 

After the soft function for the SW algorithm has been calculated, our result for the jet function can be used 
to derive a jet algorithm dependent factorized cross section for e + e~ into dijets as a function of the jet invariant mass s. 
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Note added: Right after this paper first appeared, Ref. [25( by Ellis et al. came out with a calculation of the 
full jet function. They pointed out that the treatment of the zero-bin was different in our two papers. In Appendix 
|A| we derive a systematic method to apply zero-bin scaling to phase space constraints. With this method we find the 
same zero-bin subtractions as they do. After a private communication to sort out a typographical error in the theta 
functions in their original manuscript, the results for the full jet function agree between the two papers. 
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Appendix A: Zero-bin with a Jet Algorithm 

The purpose of the zero-bin subtraction is to remove double counting between the jet and soft functions. As 
explained in Ref. [28|], there is some freedom in how to define the subtraction. This can be compared to freedom 
in choosing from different renormalization schemes, all of which remove the UV divergences but can differ by finite 
constants. Similarly, all zero-bin subtraction schemes must remove the IR divergences in a universal manner but 
may include different constants in the result. The authors of Ref. [28| advocate for a minimal approach analogous 
to minimal subtraction for renormalization. We show how their approach can be extended in a consistent way to 
apply also in the presence of phase space restrictions, such as jet algorithms. We expect there to be other consistent 
zero-bin subtraction schemes but we leave their exploration to future work. 

When dealing with phase space, it is convenient to think about the zero-bin scaling in terms of final state momenta 
instead of loop momenta. The two approaches are equivalent but the former is conceptually simpler to apply to jet 
algorithms. To begin with, we write down the Feynman diagrams for the desired process where an initial parton goes 
into a final state of several particles. The momenta of all the external particles and all the internal propagators are 
considered independent and momentum conservation at the vertices is implemented by explicit delta functions and 
integrals over the internal momenta. At any order in a s , the zero-bin contributions for a Feynman diagram are found 
by taking one or more of the collinear final state momenta pi to scale as pi ~ QX 2 . This only affects expressions where 
two different momenta are added or subtracted. In the approach advocated here such comparisons only take place at 
the vertices. Performing the integrals over the internal momenta conveys the information about the zero-bin scaling 
contained in the vertices into the rest of the diagram. The same approach can also be taken when thinking about the 
zero-bin scaling in loop diagrams that give the cross section via the optical theorem. As long as the momentum for 
every propagator is considered independent, it is possible to apply the zero-bin scaling to a single propagator and let 
the momentum-conserving delta functions convey the information about the scaling to the rest of the propagators. 
Hence, thinking about the zero-bin this way unifies the treatment of phase space integrals and loop integrals. 

In order to discuss the influence of phase space restrictions, let us consider a cross section that is fully differential 
in N variables, i.e. no integrals over physical degrees of freedom have been performed. In this case, no zero-bin 
subtraction is needed since we have full control over the momenta of all particles; we can tell which region of phase 
space they are in and hence whether they are soft or collinear. Next, let us perform some integrals that can cause the 
integrand to be evaluated in the zero-bin region while keeping the cross section differential in M variables (For 
example, think of bk as the invariant masses of the jets.) In this case, the integration region becomes a hypersurface 
in phase space, represented schematically in Fig. [5] by the blue parallelogram. The red dashed sphere represents the 
zero-bin region and the blue oval is the intersection of the zero-bin region with the hypersurface on which we have 
restricted the final state particles. The region inside the blue oval has to be removed by the zero-bin subtraction but 
there is no need to perform a subtraction outside the hypersurface. This is another way of saying that the zero-bin 
scaling does not act directly on the phase space constraints that specify the values for the bk- However, the constraints 
can be modified indirectly by the momentum-conserving delta functions. 

To state the argument more mathematically, consider two ways to perform the phase space integrals over the 
hypersurface. Either we integrate over all N degrees of freedom and use delta functions to enforce the M constraints 
or alternatively, we use the constraints to find a set of N — M independent coordinates qi and integrate over them. 
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FIG. 5: Schematic representation of the phase space. The blue parallelogram represents the hypersurface defined by phase space 
constraints, the red dashed sphere represents the zero-bin region, and the blue oval shows the intersection of the hypersurface 
with the zero-bin region. 



The two approaches can be written as 

„ N M „ N — AI 

/n d ^n^-^)= / n ( ai ) 

J i=l fe=l f=l 

If we use the right hand side of Eq. (jAlj) to evaluate the cross section, we can see that the zero-bin scaling will be 
applied to the integrand but will not affect the phase space constraints 5(bk ~ bk) directly. Similarly, we can apply 
the phase space constraints required by a jet algorithm so that the integration is performed over a region of phase 
space. Again, the zero-bin subtraction has to be performed on phase space region specified by the algorithm and 
the algorithm constraints are only modified by the effect of the zero-bin scaling on the momentum-conserving delta 
functions. 

To give a concrete example, we look at the zero-bin subtraction for the diagrams in Fig. [TJ The correspondence 
between the original momentum-conserving delta functions and the associated zero-bin scaled versions is 

5{Q-Pq ~Pg) -^5{Q-Pq) 

5 d - 2 ft - P f) s d - 2 (pf) 

S (s/Q - p+ - p+) > 6 {s/Q - p+ - p+) . (A2) 
Applying these delta functions to the on-shell condition for the quark gives 

SiPqPt + W) 2 ) 5(Q P +). (A3) 

Combining Eqs. (|A2[) and (|A3I) with the gluon on-shell condition leads to Eq. (flU|) . which together with Eq. ([12")) gives 
Eq. (US]). 

It should be cautioned that performing the zero-bin subtraction as explained above does not give the same result as 
applying the scaling p g ~ QX 2 directly to the jet algorithm constraints. In the latter approach, the gluon angle with 
respect to the jet axis would scale as A and the theta function would never be satisfied according to power counting 
since S ~ A. 

In Ref. , Ellis et al. point out that the difference between the results in our two papers could come from a different 
treatment of the zero-bin and from the different power counting that we use for S. Our calculation corresponds to 
their measured jet function with angularity tq — s/Q 2 and cone half-angle R = 5. They use power counting tq/R ~ A 2 
whereas in the present work we take tq ~ R ~ A 2 . The power counting of momenta is the same in both works. We 
carried out the investigations reported in this appendix in order to clarify these differences. Following the approach 
advocated above, we derive the same zero-bin subtractions as Ref. [25j . indicating that both calculations use minimal 
zero-bin subtractions rather than some other scheme. 

On the question of power counting of 5, if we write 5 « 2 tan 5/2, it can be seen that at leading order in power 
counting of momenta our jet algorithm constraints in Eqs. (|12p and (|13[) are equivalent to their Eqs. (5.1) and (5.3). 
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However, the terms appearing in their theta functions in Eq. (5.1) are of different order in A, whereas our theta 
functions constraining the naive contribution in Eq. (|12[) are homogeneous in power counting. For the phase space 
constraints in the zero-bin region, the scaling will not be homogeneous because, as argued above, we do not need to 
apply the zero-bin scaling to the constraints. Because Ellis et al. do not expand in A in the theta functions for the 
naive contribution, our algorithm constraints end up being identical to theirs and our Eqs. ([9J and (1111) are identical 
to their Eqs. (5.7) and (5.8) with a cone algorithm. 

The remai ning difference between the two papers in the finite part of the jet function is due to a typographical 
error in Ref. [25j where some of the theta functions have been misplaced. 2 With this fix in the theta functions, our 
two papers give the same result for the jet function. 
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